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On Minimal Spectral Factors with Zeroes and 
Poles lying on Prescribed Regions 

Giacomo Baggio, Augusto Ferrante 


Abstract 

In this paper, we consider a general discrete-time spectral factorization problem for rational matrix¬ 
valued functions. We build on a recent result establishing existence of a spectral factor whose zeroes 
and poles lie in any pair of prescribed regions of the complex plane featuring a geometry compatible 
with symplectic symmetry. In this general setting, uniqueness of the spectral factor is not guaranteed. It 
was, however, conjectured that if we further impose stochastic minimality, uniqueness can be recovered. 
The main result of his paper is a proof of this conjecture. 


I. Introduction and problem definition 

Spectral factorization is a crucial problem of many areas of systems and control theory, from 
LQ optimal control theory [8] to filtering and estimation theory [5], [6], to cite but a few. A 
seminal paper in spectral factorization theory is due to Youla [9], In that paper, a constructive 
procedure was established to compute a stochastically minimal outer spectral factor of a given 
rational spectral density <F(s) without requiring any additional system-theoretic assumption. In 
fact, Youla’s only assumptions — that are clearly necessary for the existence of a spectral 
factor — are that the real rational spectrum <f>(s) is positive semi-definite in the points of the 
imaginary axis where it is finite and features the Hamiltonian paraconjugate symmetry, i.e. 
<f>(s) = [<f>(—s)] T . Under these assumptions, Youla has established existence, and provided a 
procedure for construction, of a rational matrix function W(s ) — a spectral factor — analytic 
with its right-inverse in the open right half complex plane, such that $(s) = [W(—s)] T !U(s). 
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In [2] a discrete-time counterpart of the Youla’s result is established as a corollary of a much 
more general result that allows for the selection of the analyticity regions of the spectral factor 
and of its right-inverse. Remarkably, the latter feature is of key interest in stochastic realization 
and a-causal estimation theory, see, e.g., [3], [7]. The main result of [2], that may be viewed as 
the starting point for this note, may be described as follows. Let <h(^) be a real, rational, matrix¬ 
valued function. Assume that <3>(z) is positive semi-definite in the points of the unit circle where 
it is finite and that it features the Symplectic paraconjugate symmetry , i.e. = [ ( I>(T/;j)] t . Let 
g/ p and g/ z be regions of the extended complex plane C, compatible with symplectic structure 
(i.e. for each z G C with \z\ f 1, exactly one element of the pair (z,z _1 ) is in g/ p and 
z ^ gf p if \z\ — 1; and the same holds for g/ z ). Then, there exists a rational matrix function 
W(z) — a spectral factor — analytic in g// p and with right-inverse analytic in sf z , such that 
< L(^) = \W{1 /z)YW[z). In the case when gf p and g/ z coincide with the subset of C of the z 
that are outside the closed unit disc, we get the discrete-time counterpart of the result of Youla. 
There are, however, many other interesting situations. For example in backward filtering, g/ p is 
fixed by the system’s dynamics while g/ z is the open unit disc. 

From this general result, a very interesting question arises. In fact, when sf p = sf z = srf, it 
is not difficult to see that the corresponding spectral factor is essentially unique (i.e. unique up 
to multiplication on the left side by a constant orthogonal matrix): the key idea is that, starting 
from a reference spectral factor W(z), a second spectral factor Wi(z) must be of the form 
Wi{z) = Q(z)W(z) with Q(z) being all-pass so that if Q(z) has a pole in p, it necessarily 
has a zero in 1/p; therefore, for any non-constant Q(z), either IF, (z) or its right-inverse is no 
longer analytic in sf. On the contrary, when f we can easily obtain a spectral factor 
Wi (z) with the prescribed analyticity properties by selecting an all-pass function Q(z) featuring 
poles in C \&f p and zeroes in C \ srf z . Thus, there appears to be an inherent ambiguity in the 
choice of the spectral factor in this general case. In this paper, we show that this is in not 
the case if we further impose that the spectral factor has minimal complexity as measured by 
its McMillan degree. In fact, we will show that, under this assumption, for any choice of the 
analyticity regions gf p and sf z , the spectral factor is essentially unique. In the scalar case, this 
result is straightforward. In the general matricial case, however, a rational function can feature 
a pole and a zero in the same point so that the result appears to be quite difficult to derive 
and was left open as a conjecture in [2]. Our proof makes use of a very elegant and profound 
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parametrization of rational all-pass functions established by Alpay and Gohberg in [1], 

Paper structure. The paper is organized as follows: In §11, we review some preliminary notions 
of rational matrix theory and we introduce some ancillary results. In §111, we present our main 
theorem. Finally, in §IV, we draw some concluding remarks and we list a number of possible 
future research directions. 

Notation. In what follows, we write G T , G*, G _1 , G~ R for the transpose, Hermitian conjugate, 
inverse and right inverse of matrix G, respectively. As usual, /„ is the n x n identity matrix 
and diag[«i,..., a n ] stands for the matrix whose diagonal entries are ai,, a n . We let T : = 
{ z G C : \z\ = 1} and we denote by C := C U {oo} the extended complex plane. We denote 
by ¥[z] mxn and ¥(z) mxn the set of m x n polynomial and rational matrices with coefficients 
in the field F (we consider the two cases F = R, C). A polynomial matrix G(z) £ F[^] mxn is 
said to be unimodular if it possesses a polynomial inverse (either left, right or both). Notably, 
a square polynomial matrix G{z ) is unimodular if and only if det G{z) is a constant. Given a 
rational matrix G(z) £ C (z) mxn , we let G*(z) := [G(l/T)]*, where z is the complex conjugate 
of z £ C, and we denote by rk(G) the normal rank of G(z), i.e., the rank almost everywhere in 
z G C of G(z). The rational matrix G(z) is said to be analytic in a region of the complex plane 
if all its entries are analytic in this region. If rk(G) = m, then G~ R (z) denotes a “minimal” 
right inverse of G(z), i.e., a right inverse of G(z) whose poles coincide with the zeroes of G(z). 
Other standard notation and terminology is taken from [2]. 


II. Preliminary results 

If G(z) £ C(z) mxn , rk(G) — r < min {m,n}, by the Smith-McMillan Theorem [4, Ch.6, §5], 
there exist unimodular matrices U(z) £ C[z] mxr and V(z) £ C[z] rxn such that 


D{z) : = U{z)G{z)V{z) 

= diagl" £l (*) £ £r ^ 


( 1 ) 


where £i(z), £ 2 (^), • • •, £ r (z), fpi(z), ip 2 (z), ..., ip r (z) £ C [z\ are monic polynomials satisfy¬ 
ing the conditions: (i) £i(z) and 'ipi(z) are relatively prime, i = 1,2,... , r, (ii) £i(z) \ £i+i(z) 
and ip i+ i(z) | ipi(z), i = 1, 2,..., r - l. 1 


1 If p(z) and q(z) are two polynomials in R[«], p(z) \ q(z ) means that p{z) divides q(z) 
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The rational matrix D(z) in (1) is known as the Smith-McMillan (SM, for short) canonical 
form of G(z). The (finite) zeroes of G(z) coincide with the zeroes of £ r (z) and the (finite) 
poles of G(z) with the zeroes of G\ (z). The degree of a pole and zero at a G C (denoted by 
d p (G: a) and 8 Z (G; a), respectively) is equal to the sum of the degrees of the zero at a of all 
the i/)i(z) and of all the c jz), respectively. 2 If G(z) has no pole (zero) at a, we let d p (G: a) = 0 
(6 Z (G] a) = 0). Furthermore, if pi,... ,ph are the distinct poles (the pole at infinity included) of 
G(z), the McMillan degree of G(z) is defined as [4, Ch.6, §5] 

h 

Sm(G) := S p (G;pi). (2) 

i— 1 

Definition 1. Let G(z) G C (z) mxn , H(z ) G C (z) nxp and a G C. We say that in the product 
G(z)H(z) there is: 

1) a pole cancellation at a if 5 P (GH\ a) < 5 P (G; a) + S P (H ; a); 

2) a zero cancellation at a if 5 Z (GH; a) < S Z (G; a) + 5 Z (H ; a); 

3) a zero-pole cancellation at a if both conditions 1) and 2) are met. 

Remark 1. If rk(G) = vk(II) = n then a zero or pole cancellation at a in the product G(z)II(z) 
always corresponds to a zero-pole cancellation at a. A proof of this fact is postponed to the end 
of this section (Lemma 2). 

However, in general, 1) and 2) are not equivalent. Indeed, consider for instance the product 

2z+3 -| 

G(z)H(z) = 1 -1 (* +1 )( z + 2 ) =- 

w w L J _±_ z+1 

z+2 

and observe that there is a pole cancellation at —2 which does not correspond to a zero-pole 
cancellation at —2. 

Two special classes of rational matrices are the following ones. 

Definition 2 (Paraconjugate-Hermitian). A rational matrix G(z) G C(z) nxn is said to be paraconjugate- 
Hermitian if G(z) = G*(z). 

2 If a = oo, then we can consider the transformation z i-> A' 1 and the definition still applies by considering the degree of 
the pole/zero at A = 0 of G( A). 


September 28, 2015 


DRAFT 



5 


Definition 3 (Paraconjugate-unitary or All-Pass). A rational matrix G(z) G C(z) rxr is said to 
be paraconjugate-unitary or all-pass if 

G*(z)G(z) = G(z)G*(z) =I n . 

Remark 2. A real paraconjugate-Hermitian (paraconjugate-unitary) matrix G(z) G R(z) rxr 
is said to be para-Hermitian (para-unitary, respectively). In addition, it is worth noting that 
a paraconjugate-Hermitian matrix is Hermitian in the ordinary sense upon the unit circle, while 
a paraconjugate-unitary matrix is unitary in the ordinary sense upon the unit circle. 

A useful characterization of the class of paraconjugate-unitary matrices is provided by the 
following Lemma. 

Lemma 1. Let V(z) G C (z) rxr , Sm(V) = n, and let {aj}™ =1 be the poles of V(z) counted with 
multiplicity, then V(z) is paraconjugate-unitary if and only if it can be written as 

V(z) = UU 1 (z)U 2 (z) ■ ■ - U n (z), (3) 

with U G C rxr being constant unitary and 3 

Ui(z) := I r + (- l) p u «,gC\T, (4) 

\z-ai J 

with Pi G C rxr being an orthogonal rank-one projection. Moreover, the product in the right-hand 
side of (3) is minimal, i.e., S M (V) = S M (Ui) 4-h S M (U n ). 

Proof: See [1, Thm. 3.12]. ■ 

Remark 3. Given any decomposition of a paraconjugate-unitary matrix V(z) of the form in (4), 
we have that: 

3 We adopt the convention 1 y QL * z =: z if a, = oo. 

Z Oil 
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1) Every factor Ui(z) in (4) is paraconjugate-unitary. Indeed, by direct computation: 

U;(z)Ui(z) = I r +[ 1 ~ aiZ _ 1 - 1 ) Pi + (- 1 )P t + 


z 1 - Oli 


z — as 


-lip. 


Z 1 - Oli 


z — aj 


1 — OtiZ 1 „ 1 — Ot;Z 


= Ir- m + + — —Pi + 2 Pi - 1 _* iZ _ A - 

z 1 — a* z — a* z 1 — ai z — cti 

= L. 


2) Every pole at a t of V(z) of degree d t is accompanied by a zero of V(z) at 1 /ai of the 
same degree. In particular, if cc* ^ oo, the SM canonical form of U r {z) in (4) is given by 

1 


diag 


z — Oli 


l,z~ l/oii 


3) Since the decomposition is minimal and = 1, i — 1,..., n, it follows that 

n 

8 M {V) = Y J ^M{U i ) = n. 


i= 1 


4) Since the orthogonal rank-one projection P t in (4) can be written as Pi = v,v* with v l e 
s.t. \\vi\\ 2 = v*Vi = 1, it holds 4 


det Ui(z) = det 
= det 


Ir+ | l - sal — 1 | 


I r + 


z — ai 


1 — a.jZ 


z — ai 


~ 1 I ViV* 


= det 

1 — OiiZ 


1 i * I 1 i | 

1 + V, |-1 I Vi 


z — ai 


z — ai 

For the sake of completeness, we state and prove below two additional instrumental Lemmata. 


Lemma 2. Let G(z) G C (z) nxr and H(z) G C (z) rxm with rk(G) = rk (H) = r. If G(z)H(z) 

has a zero or pole cancellation at a G C, then G(z)H(z) has a zero-pole cancellation at a. 

Proof: Assume that G{z)H{z) has a pole cancellation at a G C (the proof for the case of 
a zero cancellation ataGC goes along the same lines). 


4 In the derivation we exploit the fact that if A € C nxr and B G C rxn , then det(/„ + AB) = det(/ r + BA). 
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Let D(z), D'(z ) G C (z) rxr be the SM canonical form of G(z), H(z), respectively. We can 
write 

G(z) = C(z)D{z)F{z ) and H(z) = C\z)D\z)F\z) 


with F(z), C'(z ) G C[^] rxr , C(z) G C[z] nxr and F\z) G C [z] rxm unimodular matrices. Hence 
the product G(z)H(z) can be written as 

G(z)H(z) = C(z)D(z)F(z)C\z)D'(z)F'(z) 

where M(z) := F(z)C'(z) G C[^] rxr is unimodular. Notice that, by virtue of the unimodularity 
of C{z) and F'(z), the SM canonical form of G(z)H(z), denoted by A(z), coincides with that 
of D(z)M(z)D'(z) (see [4, Ex.6.5-6]). Moreover observe that, since rk(G) = rk(if) = r, then 
rk (DMD') = r. Therefore, by taking determinants, we have 


detA(z) — cdet D(z) det D\z) 

n(z) (z — Q;) <5 dG';aO+<5zG H '; Q ) 


= c- 


(5) 


d(z) (z — Q;) <5 p(Aa)+i5p(^;a) 
with n(z) and d(z) relatively prime polynomials s.t. n(a) ^ 0, d(a) ^ 0, and c G C, 0. On 
the other hand, since A ( 2 ) is the SM canonical form of G(z)H(z), we get 

n(z) (z — a ) Sz ( GH;a ) 


det A(^) = c 

Hence, a comparison of (5) and (6) yields 


d(z) (z — ay P ( GH ; a ) 


( 6 ) 


5 P (GH- a) - 5 P (G ; a) - 5 P (H ; a) = S Z (GH ; a) - S Z {G\ a) - 5 Z (H ; a). (7) 

Since, by assumption, G{z)H{z) has a pole cancellation at a, the left-hand side of (7) is strictly 
negative. This in turn implies that the right-hand side of (7) is strictly negative, i.e. G(z)H(z) 
has a zero cancellation at a. From this fact the thesis follows. ■ 


Lemma 3. Let G(z) G C(z) nxr and H(z) G C(z) rxm with rk(G) = rk (H) = r. If G(z) and 
H(z) have no zeroes at o. G C then 

6 P (GH; a) = S P (G; a) + 5 P (H; a). 
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Proof: By following verbatim the first part of the proof of Lemma 2, we arrive at the 
expression 

det A(z) = det D(z) det M(z) det D'(z ) 

= cdet D(z) det D\z), c £ C, c f 0. (8) 

Since by assumption G(z) and II (z) have no zero at a, then D(z) and I)'(z) have no zero at 
a. Furthermore, A(z) has no zero at a. This fact can be seen by taking the inverse of A(z), 
namely 

A _1 (A) = D r ~ 1 (z)M~ 1 (z)D~ 1 (z), 

and by noting that the latter has no pole at a, since the entries of D^(z), M~ l (z) and D , ~ 1 (z) 
do not have any pole at a. This in turn implies that S p (D; a), d p (I)': a) and 5 p (A; a) coincide 
with the degree of the pole at a in detL>(A), det D'(z) and det A(^), respectively. Hence, by 
summing up all the previous considerations, we get 

5 p (GH;a) = 5 p (A;a) 

= <5 p (det A; a) 

= (5 p (det D\ a) + 5 p (det D'\ a) 

= 8 P (D-, ot) + 8 P (D'-, a) 

— 8p(G ; o) T 8 P (H; cr) 

which concludes the proof. ■ 

Remark 4. Notice that Lemmata 2 and 3 still hold when a = oo. As a matter of fact, in this case, 
we can apply the change of variable 2 A^ 1 and then consider the (degree of the) zero/pole at 

A = 0 in G( A) and H( A). 


III. The main theorem 

Before proceeding with the proof of the main Theorem, we introduce some preliminary 
definitions. 

Definition 4 (Spectrum). A para-Hermitian rational matrix <f>(A) e M(^) nxn is said to be a 
spectrum if ^(e- 7 ^) is positive semi-definite for all 00 £ [0, 2n) such that <f>(e- 7 ^) is defined. 
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Definition 5 ((Stochastically minimal) Spectral factor). Given a spectrum <3>(z) e R(z) nxn , a 
real matrix-valued function W(z) satisfying 

$(*) = W*(z)W(z), 

is called a spectral factor of $(z). Moreover, the spectral factor W(z) is said to be stochastically 
minimal if 

Sm(W) = 

Stochastically minimal spectral factors correspond to solutions of minimal complexity (i.e. 
minimal McMillan degree) of the spectral factorization problem. Indeed, a spectral factor W(z) 
of ^{z) cannot have McMillan degree smaller than \ m*)- 

Definition 6 ((Weakly) Unmixed-symplectic). A set srf C C is unmixed-symplectic if 

srf U srf* = C\ {z£C : \z\ = 1}, and srf D sf* = 0, 

where sf* = {z : z _1 £ srf}. The set sf C C is weakly unmixed-symplectic if 

srf U = C, and srf D sf* = {z £ C : \z\ — 1}, 

The following Theorem is the main result of this note. 

Theorem 1. Let <f>(^) e M(z) nxn be a spectrum with rk(<3>) — r < n, r f 0. Let W(z), Wi(z) £ 
R(z) rxn be such that 

1) W(z) and W\{z) are spectral factors of <&(;?), i.e. <f>(^) = W*(z)W(z) = W*(z)W\(z)\ 

2) W(z),Wi(z) are analytic in and IL R (z),Wf R (z) are analytic in stf z , where sf v , 
are weakly unmixed-symplectic regions; 

3) W(z) and W\(z) are stochastically minimal, i.e. S M (W) = S M (Wi) = |<5 M -(<f>). 

Then, W\(z) = TW(z) with T £ M rxr constant orthogonal. 

Proof: Before illustrating the details of the proof, we outline the key steps in order to 
provide a road-map that may help the reader. 

i) We consider the para-unitary function T(z) satisfying W\{z) = T{z)W{z) and we show 
that T(z) must have no poles and zeroes in the region (T sf z . 

ii) We then assume by contradiction that T(z) is non-constant and, more precisely, that T(z) 
possesses poles both in \ srf z and in srf z \ srf v . 
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iii) We decompose T(z) according to Lemma 1 and, by exploiting the properties of this decom¬ 

position, we show that for each pole a G sZ v \ srf z of T(z) there is a zero-pole cancellation 
both at a and at 1/a in the product T(z)W(z). Hence, we arrive at the contradiction that 
there exists a spectral factor of $(z), say W(z), such that 5m(W) < Since this is 

not possible, we conclude that T(z) must have no poles in the region srf v \ s/ z . 

iv) Finally, we exploit the fact that, by point 3), W(z) and W\ (z) are stochastically minimal 
spectral factors to conclude that T(z) must have no poles in the region srf z \srf v . This implies 
that T(z) is a constant and orthogonal matrix. 

We now describe the details. Consider the matrix 

T(z) := Wi(z)W~ R (z). 

By taking into account Property 1), it is immediate to see that T*(z)T(z) = /, i.e. that T(z) is 
para-unitary. Moreover, since the inverse of T(z) is given by 

T*(z)=T-\z) = W(z)Wr R (z), 

it follows that T(z) is analytic with its inverse in srf z n sf v . Now observe that 

W 1 (z)=T{z)W{z). (9) 

To see this, set Z(z) := W\{z) — T(z)W(z). By employing again Property 1), it is immediate 
to see that Z*(z)Z(z) = 0 so that Z(z) is identically zero in the unit circle and, eventually, 
Z(z) — 0. We need to show that T(z) is constant. 

Assume, ab absurdo, that T(z) has McMillan degree d with poles ai,...,a n of degree 
mi, ..., m n (d — mi + ■■■ + m n ), respectively, s.t. a 1 ,...,a t G sZ p \srf z and a t + 1,..., a n G 
g/ z \ srf v . In what follows we assume that 7 ^ 00 for i — 1,..., n. As a matter of fact, if this is 
not the case, we can always find a suitable Mobius transformation z (->• f(z) such that T(f(z )) 
has only finite poles. Thus, by considering this transformation, the argument in the proof still 
applies. 

By exploiting Lemma 1, we can decompose T(z) as 

T(z) = UU aitl (z) ■ ■ ■ U aijmi -i(z)U an (z) ■ ■ ■ U a2 (z)U aumi (z), (10) 
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with U E C rxr constant unitary and 

U aiJ (z) := I r + _ A P ( 11 ) 

\z-ou J 

U Ql (z) :=U Qitl (z)---U a ^ mi (z), (12) 

with i — 1,_, n. j — rn u . . m n , and P l} e C rxr being an orthogonal rank-one projection. 

Now, we can rearrange (9) in the form 

• ■ ■ K^U-W^z) = U an (z )... U„(z)U mm (z)W(z). (13) 

Notice that the left-hand side of (13) is analytic in \ srf z with (right) inverse analytic in 
■v/ z \ It follows that the right-hand side of (13) must be analytic in xy/ p \ with (right) 
inverse analytic in srf z \srf p . By rewriting the right-hand side of (13) in a more explicit way, we 
obtain 

Ua n (z) ■ • • U aa {z)U aumi {z)W{z) ( = U an (z) • • • U a2 (z) (ir- P 1>mi + 1 - — P^A W(z) 

V z-a I ) 

= u an (z) • • • U a2 (z)^—^P hmi W(z) + A(z) 

z — a i 

where A(^) := U an (z ) • • ■U a2 {z){I r — Pi }ini )W(z) has no pole at oi\. In fact, an e s# v \ srf z 
so that W(z) does not have a pole at a i. The minimality of the factorization of T(z) in (10) 
implies that 

(l-a^z p \ W{z) 

\z-a i J 

must have a zero-pole cancellation at o>\. This fact needs a detailed explanation. 

First, define U ies := U an (z ) • • • U a2 (z ) and notice that, since the factorization of T(z) in (10) 
is minimal, the matrix 

R(z) := U res (z)^—^Pi, mi 

z — a i 

has a pole at a\. In fact, a pole cancellation at a\ in R(z) would imply a pole cancellation at 
in U an (z ) • • ■U a2 (z)U ai ,m 1 {z), yielding that the degree of the pole a\ in T(z) is less than 
m\. However, this is not possible since, by Lemma 1, the factorization in (10) is minimal. Now, 
since P\, mi is an orthogonal rank-one projection, there exists a unitary matrix 0 e C rxr such 
that 

Q*Pi, mi Q = diag[l,0,... ,0], 


September 28, 2015 


DRAFT 



12 


Since Q is constant and nonsingular, also 

R(z) ■= R(z)Q = U res (z) ] —^QQ*P 1 , ini Q 

z — a i 

= U ies (z)Qdmg -——, 0,..., 0 (14) 

z — a i 

has a pole at a\. More in detail, at least one entry in the first column of R(z) possesses a pole 
at «!, while all the other columns are identically zero. Now consider A(z) := R(z)W(z) 

A{z) := R(z)W(z) = R(z)W(z ), (15) 

where W(z) := Q*W{z). As already observed, A(z) is analytic in zz/ p \zA z . Therefore, by taking 
into account that at least one entry in the first column of R(z) possesses a pole at ay, while 
all the other columns are identically zero, it is immediate that every element in the first row of 
W(z) has a zero at ay or is identically zero. Then, 

Q diag -——, 0,...,0 W{z)=(- — ^P hmj \w(z) 

_ z — a i J \ z — a \ ) 

has no pole at a^. This implies that also U aitTni (z)W(z) has no pole at oil so that in the 
product U aitrni (z)W(z) there is a pole cancellation at a,\. Eventually, since U nurrl] (z) has full 
(column-)rank and if" (z) has full row-rank, by Lemma 2, we can conclude that in the product 
^ f a 1 ,m 1 (z)W(z) there is a zero-pole cancellation at a\. 

By replacing (13) with 

Wi a (z)VV n (z) ■ ■ ■ W = W- R (z)U’ umi (z)U’ m (z) ■ ■ ■ Uljz), 

we can repeat almost verbatim the previous argument in order to conclude that W~ R (z)U* l (z) 
must have a zero-pole cancellation at 1/ai, or, equivalently, U aijmi (z)W(z) must have a zero- 
pole cancellation at 1 /ay. 

The zero-pole cancellations at ot\ and at 1/ay in the product U aijTni (z)W(z) imply that 

S M (U ai , mi W) < S M (W) = 

Indeed, let p u ... ,p h e C be the poles of W(z) s.t. 7 ^ 1/ai for all % = 1,..., h. Since 
Ucn,mi(z) is analytic together with its inverse in C \ {ai, 1/ay}, it holds 

S P (W] Pi )=S p (U aitmi W] Pi ) 
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for all i — 1,..., h. Moreover, by the zero-pole cancellations: (i) (z)W(z) has no pole at 

ai, and (ii) S p (U Qlyrni W; 1/ai) < 5 p {W\ 1/ah). Therefore 

h 

S M (U ai , mi W) = S P (W;+ S p (U ai:mi W] l/a\) 

1=1 

h 

< J2 s p( w '’Pi) + s p ( w ; V«i) = 

i=l 

But this is clearly not possible since, by point 3), W(z ) is a stochastically minimal spectral 
factor. Therefore, U ai:7ni (z) must be a constant unitary matrix. 

The previous reasoning still applies for all the other factors of T(z) having a pole at a t , 
i = 1,..., t, yielding that rri t — 0 for all i = 1,,.., t, i.e., T(z) has no poles at i — 1,..., t. 
It remains to show that T(z ) has no pole at a t+ 1 ,..., a n . To this aim, we have 

Wi(z) = T(z)W(z) 

and since all the poles of T(z) lie in srf z \ by Lemma 3, we have 5 p {Wi\ «*) = 5 P (T\ ctj) + 
5 P (W] ttj) for all i = t + 1,..., n, while for all the other poles p. t , i = 1,..., h, of W(z), 
8 p (I'Li; Pi ) = 8 p (W-,Pi). This implies that 

n n h 

S M (W i )= ^ 6 p (T-Oi)+ Y, 5p{W- ai ) + Y,^ W ^^ 

i=t -\-1 i=t- 1-1 i=l 

n h 

> ^ 5p{W ] a i ) + Y J ^v(W ] p i ) 

i=t+1 i=1 

= s M (w), 

which, by virtue of the stochastic minimality of Wi(z), leads to a contradiction. Hence T(z) 
must have no poles at a t + 1 ,..., a n . 

To conclude, we have shown that T(z) has no poles and hence no zeroes, due to the fact 
that T(z) is a para-unitary matrix. Therefore, since it has real entries, T(z) must be a constant 
orthogonal matrix. ■ 


IV. Conclusions 

In this paper we have analyzed uniqueness of the solution of spectral factorization problem 
with prescribed dynamical features. If we restrict attention to solutions of minimal complexity, the 


September 28, 2015 


DRAFT 



14 


solution is indeed essentially unique. The proof is based on the parametrization of discrete-time 
all-pass functions provided in [1] and on some preliminary results on rational matrix functions 
that we have established and that may be of independent interest. 

We believe that similar techniques may be employed to derive the continuous-time counterpart 
of this result. Indeed, in [1] a parametrization of continuous-time all-pass functions is also 
provided and the rest of the procedure appears to be adaptable to the continuous-time case as 
well. 
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